Decay of Langmuir wave in dense plasmas and warm dense matter 
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The decays of the Langmuir waves in dense plasmas are computed using the dielectric function 
theory recently proposed [lj. Four cases are considered: a classical plasma, a Maxwellian plasma, 
a degenerate quantum plasma, and a partially degenerate plasma. The results obtained suggest 
considerable deviations from the conventional Landau damping theory. Its implication on the x-ray 
Raman compression in dense plasmas or warm dense matter is discussed. 
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I. INTRODUCTION 

The decay of the Langmuir wave is one of the most im- 
portant physics in plasmas related to many applications 
[l]-[7jj. While the Landau damping theory is accurate to 
predict the plasmon decay in ideal classical plasmas, it 
is inadequate in dense plasmas and warm dense matter. 
For an example, in metals, the experimental measure- 
ment shows that the decay of the long- wave length plas- 
mon is higher than that predicted by the Landau damp- 
ing theory Q. Previously, a new theory ||, which is a 
natural extension from the one in the condensed mat- 
ter [I3-[l2j], has been proposed to compute the plasmon 
decay in dense plasmas. 

In this paper, we provide detailed calculations of this 
theory [9(. We consider four cases: a classical plasma, a 
Maxwellian plasma, a degenerate quantum plasma, and 
a partially degenerate plasma. The decays in the limit, 
where the wave vector goes to zero, are finite in all cases. 
The decay rate is lower in finite-temperature plasmas 
than degenerate plasmas. We then discuss the implica- 
tion of the result on the x-ray Raman compression where 
dense plasmas or warm dense matter are used as a com- 
pressing media. 

This paper is organized as follows. In section II, a 
new theory [i| is introduced to predict the plasmon de- 
cay in dense plasmas and warm dense matter. In section 
III, we compute the decay rate implied by the theory 
Q in the limit when h = 0, where h is the Frank con- 
stant. In section IV, we compute the plasmon decay for 
a Maxwellian plasma when ft ^ 0. In section V, the plas- 
mon decay is computed for the quantum electron plasma 
(h 7^ 0) when T e — 0. In section VI, the plasmon decay 
is computed when T e ^ and H ^ 0. In section VII, 



the summary is given, and we discuss the implication of 
the result obtained on the x-ray Raman compression in 
dense plasmas. 



II. DIELECTRIC FUNCTION THEORY 

Previously, it is shown that the dielectric function in 
dense plasmas under a wave with wavevector k and an- 
gular frequency lu has additional term [9|: 



47re 2 47re 2 
e(k, w) = 1 + -p-a rpa (k, u) + -p-a do nsc(k, u), (1) 

where adonso is the special term that is significant in dense 
plasmas. The a rpa is the well-known Lindhard suscepti- 
bility 



ri 3 q /(|k + q|)-/(q) 
(2tt) 3 Hu - E(k + q) + E(k) ' 



(2) 



where / is the occupation number and E(q) = H 2 q 2 /2r 



The adenso is computed [9( as 

d 3 s , s f d 3 q U 2 (s) 

x3»W) 



"dense (k,w) 



(27r) 3 ""^ 7 7 (2tt) 3 A 2 (q,k,s) 
/(|q + k + s|)-/(q) 
hui- S(|q + k + s|) + £:(q)' 



(3) 

where U(s) = 47rZje 2 /(s 2 + kf) is the screened ion- 
electron potential with screening length l/k s , Ui(s) is 
the ion-ion structure factor rii(s) — (J"^- ■ exp(is • (Xj — 
Xj))/F, A(q, k, s) is defined as 

A^fq.k.s) = — — - — 

^' > E(q)-E(\ q + s\) 
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£(|q + k|)-£(|q + k + s|) 



(4) 
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The imaginary part of a<iense * s given as 



f d 3 s f d 3 q lP{s) 

Im[a dcnsc ] = J ^3«i(8) y (27r)3i42(q>k> 



d 3 q t/ 2 (s) 



X7T(5(^ - £(|q + k + s|) + E(q)) 
x(/(|q + k + s|)-/(q)). 



(5) 



III. CLASSICAL PLASMAS: h = 

In the limit fi. = 0, the A in Eq.Q is simplified as 

k ■ s , 



A- 1 (k,s,q 2 ) 



(6) 



where 



^(JfejSjCOS^)) = 1 



fc 2 + s/ccos(#)\ 1 



|s + k| ; 



and cos(0) is the angle between k and s. From Eqs. ([3]) 
and ((6|), the Eq.© can be simplified as 



Aire 



■Im [adonsc] = J 



n l (s)d 3 s . , _ , 2 

/ 2 A 3 (K + (k,s) K (k,s)j 



ks 



47rZ.,e 2 



W x2a/( w /|k + s|) 



where fcdc = 47rn e e 2 /T e is the Debye screening length. 
By using the spherical coordinate, it can be re-casted as 



47re" 



Im [ttdense] 



V2tt f Z, 



where the logarithmic factor, log(A), is given as 



log(A) 



ds 



+i 



.2 ,M S ) 



u du- 



o |k + s|y_! 

x (k~k + ) 2 exp(— 



s 2 + fc 2 



2|k + s|V r 



0. 



and u = cos(#), w tc = \/T e /m e and is the ion average 
density. In particular, in the limit k — 0, = 1, the 
logarithm factor is given as 



i / a \ ds2 
log(A) = / — - 
Jo s3 



ds 2 rii(s) 



exp(- 



This integration is divergent in the s-integration due to 
the lack of the high cutoff. We introduce the high cut- 
off as Shigh = l/&c where b c = e 2 /T e is the closest ap- 
proaches. Since the decay rate is proportional to n. 
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FIG. 1: Damping rates of classical 400 eV hydrogen plasma 
(Classical), degenerate plasma (Degenerate), and partially 
degenerate plasma (Partially-degenerate), where 10 24 /cc < 
n e < 10 26 /cc and k S> &f- The electron temperature of the 
partially degenerate plasma is 0.5Ef(n e ). The damping rate 
7 is normalized by u; pc . 



and T e , its rough interpretation should be the ion- 
electron collision. In Fig. [TJ we show the decay rate, in 
the limit k = 0, normalized by CiJpQ clS cl function of the 
electron density for the long wave length plasmon, where 

CJ pc = (47T71 e e 2 /TO e ) 1 / 2 . 



IV. MAXWELLIAN PLASMA: h 

In this section, we consider the case when the electron 
distribution, /, is given as a Maxwellian and ignore the 
Fermi degeneracy. In this case, the susceptibility of the 
Lindhard function in Eq.Q is given as 



h(io — k • v) 



(8) 



where f m = I/(2vru t 2 c ) 3 / 2 exp(-(v 2 /2u te )), and v = 
hc\/m e is the electron velocity. In case when k is not 
zero, the A in Eq.© is dependent on q. The rather 
complicated expression for A is given as 



h 2 k- s 



1 



Hk- (q + s/2) 



m e {huj) 2 

hk - (q-s/2) 



mui 
l 



A- 1 (k,s,q) 



A cannot be taken out of the q-integration in Eq. ([3]) due 
to its dependency on q. For a given s and k, we do the q- 
integration first and then s-integration later. For simplic- 
ity, assume k = (k x ,0, 0) and s = (s cos(#), s sin(0), 0). 
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Let us represent the q = q±xi + q2%2 + 93^3, where 
i\ = (fc + scos(0),ssin(0),O)/|fc + s|, X\ = (— ssin(0), 1 + 
cos(0),O)/|fc + s\, and £3 = (0,0,1). In this coordinate 
system, A is only dependent on fc, s, 9, and q2'. 



A 1 (k,s,q 2 ) = 



k- s 



(9) 



where 



fc 2 + skcos(9) 

WTW 2 



h sk sin(0)g2 _j_ sk cos(0) , 
fn e uj |s + k| 2 



The real part of a in Eq.© is still very complicated, but 
the imaginary part given in Eq.© could be simplified 
due to the fact that the delta function will eliminate qi 
integration and the integrand of q^ can be done easily: 



d 3 g U 2 {s) 

(27T)3^(q J k,8) 



Im 



/(q + k + s)-/(q) 



hus - E(q + k + s) + £(q) 




FIG. 2: Damping rate of 36 eV hydrogen plasma as a function 
of q, where n e — 10 24 /cc. The x-axis is the wave vector k/k,F 
and the y-axis is damping rate ^y/uj pc . 



(k-s) 2 U 2 (s) 



dq 2 (f m {Q + ,q2,uj) - f m (Q ,g 2 ,w)) 



To begin with, the susceptibility a rpa of the free electron 
plasma has been computed by Lindhard [l3| and is given 
as 



ft 2 |s + k 



(10) 



a r p a (k,w) 



3n„ 



m e vi 



-h(z,u), 



where Q ± (\k + s\) = m e uj/h\s + k\ ± \k + s\/2. Note that 
f m is the Maxwellian distribution with q% having been 
integrated out: 



/m(<7l,<72,w) 



2?rfc 



— ex P( ^2—)' 



le 



2q 



where q tc — y/m e T e /h 2 . The right-hand side of Eq. (fT0|) 
is only a function of fc, s, 9. We could do the s-integration 
in a spherical coordinate to have the final expression as 

f 2ns 2 dfj,ds (ks^i) 2 ( AnZie 2 x 
Im[a(k,u,)] = J __„,(,)__. (^=3 

/OO 
^92 (fm(Q + ,<l2u) - / m (Q~,<?2, w)) 
"OO 

■( K +0 2 . (ii) 



ft 2 |s + k 



In Fig. [21 we plot the decay rate as a function of fc when 
n e = 10 25 /cc. 



V. DEGENERATE CASE: ft / AND T e = 

Consider a dense plasma, where electrons are com- 
pletely degenerate. The Fermi-energy is given as Ep = 
36.4 x (n/n 2 4) 2 / 3 eV where n 2 4 = 10 24 /cc, and as long as 
T e <C Ep, the computation in this section will be valid. 



where vf — \/2Ep/m e is the Fermi velocity, z — k/2kp, 
u = ui/kvp, and h = h r + i/i,-. The real part of h is given 



as 



1 1 / , kO\ , (\z — U + II 

|g + u + 1| 
|z + M+l| 



+ j- z (l-(z + u) 2 ) log( 



The imaginary part of ft, is given as 

hi — %u, if\z + u\ < 1, 

^(1 -(z-u) 2 ), if|z-u| < 1< + 



0. 



if 1 2: - u| > 1. 



In the limit where fc = 0, while computing the imagi- 
nary part of adense, A can be approximated as A^ 1 = 
(k • s)/m e uj 2 , due to the delta function in Eq.©. A be- 
comes independent to q and it can be taken out from 
the q-integration in Eq.©. Using Eq. ([2]), Eq.© can be 
simplified to 



Im[a dense (k,w)] = 



d 3 s (k-s) 2 (47rZ,e 2 ) 2 
(2tt) 3 m 2 w 4 (s 2 + k 2 tf ) 2 

xni(s)Im [a- rpa (k + s, a;)] . 
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The above equation can be further simplified to 
AE%Zi{k F e 2 ) 2 (4tt) 2 f d 3 s 



Im[a d on SO (k,w)] = 



(hw) 4 

(k-s) 2 

k 2 (s 2 + k 2 f ) 2 



3tt 2 

lm[a ipa (\k + s\,u)] ,(12) 



where we assumed rij(s) = rij. For U(s), we will use 
the screening potential J7(s) = 47rZie 2 /(s 2 + fc 2 y), where 

fc f / = ^3u> 2 e /v' F is the Thomas- Fermi screening length. 

In Fig. [TJ we plot 7/w pe as a function of electron density 
n e using the above equation (Degenerate). 

In the case of k ^ 0, the above simplification is not 
possible. Using Eq.Q, we need the q\, q 2 and q% inte- 
grations. The real part of the susceptibility in Eq. ((3|) is 
complicated but the imaginary part given in Eq. is 
simpler due to the delta function. The q\ integration is 
eliminated by the delta function and, the (^-integration 
can be done since the integrand is independent of q%. 
We will not present the detailed steps, but, after tedious 
manipulation, the q integration can be simplified to an 
1-dimensional integration given by 




FIG. 3: Damping rates of a degenerate hydrogen plasma as a 
function of the wave vector when n e = 10 26 /cc. The x-axis is 
wave vector k/kF and the y-axis is the damping rate 7/wpe. 



VI. PARTIALLY DEGENERATE CASE: h 
AND T e 



U 2 {s) 



(27r)3A 2 (q,k,s) 



Im 



/(q + k + s)-/(q) 



hw - E(q + k + s) + E(q) 



(h 2 k-s) 2 U 2 (s ) ^l-(z-u) 
ra^ihbj) 4 



kf l — (z — u) 2 



_dq 2 



47r 2 ft 2 ls + k| 



xk + k 1 (i? (k F ,z,u,q 2 ) - R(kf,z,u,q 2 )) , 



(13) 



where R ± (kp, z,u,q 2 ) — \Jk 2 F {l — (z ± u) 2 ) — q 2 , z — 
\k + s\/2kF and u = ujjkvp. The right-hand side of the 
above equation is only a function of k, s, 9. We could do 
the s-integration in spherical coordinate to have the final 
expression as 



.. 27T S 2 ^ S {h 2 k.^) 2 U 2 { S ) 

lm(a(k, uj)) = I — — — - — rii{s)- 



(27T) 



m 2 (huj) 4 



,dq 2 



-K K 



X (R~(kp,z,u,q 2 ) - R(kf,z,u,q 2 )) , 



(14) 



In Fig. [3J we plot 7/w pc in Eq. ([T4"|) as a function of k 
for a hydrogen plasma with the electron density of n e = 
10 26 /cc. 



Due to the partial degeneracy and quantum diffrac- 
tion, this case is very hard to treat. In Eq.©, the 93- 
integration cannot be simplified due to the partial de- 
generacy. However, there is some simplification achiev- 
able due to the Dharma-Wardana's technique [T3|. For 
partially degenerate electrons, the occupation number is 
given as 



f(T e ,n e ,E) 



exp 



E - /j,(n e ,T e 



(15) 



where E is the electron kinetic energy, /1 is the chemical 
energy. Dharm-Wardana showed that the / in Eq. (fT5")) 
could be expanded as follows. 



f(T e ,n e ,E) = 



2T e 1 + cosh 



fj, — S 



xg(0,n m (S) : E)dS, 



(16) 



where n m (S) = (I /3ir 2 )(2m e S / ti 2 ) 3 / 2 is the electron den- 
sity whose Fermi energy is S, and g(0, n m (S), E) is the 
occupation number for the zero temperature plasma with 
the electron density of n m (S). The Eq. (TTtll) makes it pos- 
sible to express the dielectric function of the electrons as 
the sum of the dielectric functions of the zero tempera- 
ture electrons. For an example, the susceptibility of the 
Lindhard dielectric function for electrons with a non-zero 
temperature could be expressed as an integral of the sus- 
ceptibility of the zero-temperature electrons: 
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p OG 

a Tpa (k,uj,T e ,n e ) = / a lpa (k,uj,0,n m (S)) 
Jo 



2T e (l + cosh( ;j( " e ^ ) - S )) 



dS. 



We can apply the same technique. First, define, 



B(K^O,n e )=J (27r)3A2(qjk) 



s) 



xlm 



/(q + k + s) -/(<?) 



ftw - J5(q + k + s) + E(q) 



The B(k,w,0,T e ) is given in Eq. (13]). Then, using Eqs. 
([5|), (fT3|) and ([TrJl) . the Im(a<iense) is given as 



Im(adonsc) = 
d 3 s 



o 2T e (l + cosh(^^^)) 



(27T) 



ni(«)fl(|k + s|,w, 0,7^(5)) 



dS 



(17) 



In particular, in the limit of k — or with the 
approximation that k = 1 in Eq. ([9]), the damp- 
ing rate could be obtained from Eq.(fl"2"|) by replacing 
e ipa ,(\k + s\,w,0,n e ) - 1: 



e rpa (T e ,n e ) - 1 = / (e rpa (0,n m (5)) - 1) 
Jo 



2T e (l + cosh( 



(18) 



In Fig. [1] the decay rate of the long wave-length plasmon, 
using Eq. (fl~7f , is plotted as a function of n e for a hydrogen 
plasma. 



VII. CONCLUSION 

The Landau damping theory is inadequate in dense 
plasmas; The decay rate observed from the electron stop- 
ping experiment in metals is finite. Previously, a new 
theory Q has been proposed to predict more accurately 
the plasmon decay in dense plasmas. In this paper, we 
provide the detailed calculations of the theory [9| for var- 
ious regime of dense plasmas. 

First, we consider classical plasmas where H = 0. The 
damping rate resembles the electron collision rate, which 

is proportional to rii and T e 3 ^ 2 . According to the theory, 



the 7/wp is, for a fixed temperature, an increasing func- 
tion of n e as shown in Fig. (JlJ . In the integration, there 
is a divergence in the high k, which could be avoided 
by introducing the cutoff given by the closest approach. 
The theory is be valid if the de Broglie wave length is 
smaller than the closest approach. Second, we consider 
a Maxwellian plasma. Contrary to the classical plasma, 
the high cutoff is not needed since it is provided by the 
thermal De Broglie wave length. This theory is useful 
for hot dense plamsas, where the quantum diffraction is 
not negligible, but the quantum degeneracy is. Third, 
we consider a completely degenerate plasma. Due to the 
degeneracy, 7/cj is a decreasing function of n e as shown 
in Fig. ([T]), which is an opposite case with the one in 
classical plasmas (as the electron density gets higher, the 
classical theory breaks down since it fails to take into ac- 
count the degeneracy) . The rate computed from Eq. f| L2[) 
is higher by 5 times than the experimental data [8j or the 
calculation by Sturm [Toj . This is not surprising since 
the Umklapp process in metals are very different from 
those in dense plasmas where the ion lattice structure 
is absent. Lastly, we consider the partially degenerate 
plasma. The damping is reduced in comparison to the 
degenerate case. This is natural since the ion-electron 
collisions are reduced in the partial degenerate case com- 
pared to the degenerate case [l5|, [l6[ ■ This regime is the 
most important in the application of the Raman compres- 
sion in dense plasmas or the warm dense matter. For an 
example, consider a hydrogen plasma with the electron 
density of 10 26 /cc and the temperature 400 eV. This 
case is shown in Fig. [TJ The classical theory is not valid 
in this regime and it is necessary to use Eq. (|17[) . 

The result obtained in this paper could have many im- 
plications for processes involving the Langmuir wave in 
dense plasmas, warm dense matter and metals. While 
the theory might need more refined adjustments such as 
the strong correlation, the local field correction, and the 
exchange interaction, a complete profile of the damping 
rate is, if still rough, now available as a function of the 
temperature, density and the wave vector k, which could 
be important for various processes in dense plasmas [l|- 

BiEi. 

Now, we discuss the implication of the result discussed 
above on the x-ray Raman compression. In dense plas- 
mas, the higher the temperature is, the lower the inverse 
bremsstrahlung is. It is shown previously that the Lan- 
dau damping of the plasmon generated from the ponder- 
motive potential of a pump and a seed is greatly reduced 
in dense plasmas due to the electron quantum diffrac- 
tion while the decay of the background noise plasmon 
could be very heavy due to a high electron temperature. 
This suggests that the premature pump depletion from 
the background BRS is easy to suppress while the BRS 
compression is still possible 19]. If this is the case, the 
inverse bremsstrahlung and the FRS are the most im- 
portant physics to check for the plausibility of the BRS 
compression. The optimal parameter regime would be 
determined by this consideration. 
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For T e = 0, 7/wpe is shown, in this paper, to decrease 
with an increasing temperature, and for a fixed finite 
temperature, it is shown to decrease with an increasing 
density. If the FRS is too severe for T e = 0, the FRS 
could not be contained for T e > unless the density n e 
is higher. Now, assume that the FRS is contained when 
T e — for a given n e . As the temperature increases, the 
FRS becomes stronger. Choose the maximum J"™ ax such 
that the FRS is still tolerable. This T e max might be the 
optimal temperature since the inverse bremsstrahlung is 
minimal among the parameter regime where the FRS is 
contained. 

For a high electron temperature, the FRS might be 
weak due to the enhanced Landau damping instead of 
the damping from the Umklapp processes while the BRS 
is still strong due to the reduced Landau damping from 
the band gap [l9j . In this case, the optimal temperature 
would be the maximum temperature among the param- 
eter regimes where the BRS plasmon is still a collective 
mode. Whether the FRS is contained by the Umklapp 
process or the Landau damping will be important factor 
in the determination of the optimal physical parameter 
regime, which depends on the pulse duration and inten- 
sity. We leave this question to a future researches. 



The ionic structure factor depends strongly on the ion 
temperature in the warm dense matter as the ions ex- 
perience various phase transitions with the increasing 
temperature. As shown here, the plasmon damping is 
sensitive to the ionic structure factor. This strong de- 
pendence of the plasmon damping on the ionic structure 
factor could be useful in the Raman compression, and 
it could also serve as a diagnostic in the X-ray Thom- 
son scattering in the warm dense matter [2(| [2l[ . . In 
the case of the warm dense matter, the experiment to 
measure the damping rate might be readily available us- 
ing the thin heated foil experiment ||. The effect of the 
phase transition on the plasmon decay is theoretically 
challenging, but very important for its practical appli- 
cation in the BRS x-ray compression and for other pro- 
cesses. It would be interesting to see how the damping 
changes as the ionic structure factor changes. 

In addition to the theory proposed here, the physical 
processes in dense plasmas could be very different from 
those of rare dense plasmas due to the electron degener- 
acy [H, EH 1 the electron quantum diffraction [19| and 
the band gap [13, [3- The study and applications of 
those new physics could be exciting and interesting. 
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